o 
o 

< 



^ 



DESY 06-029 ISSN 0418-9833 

April 2006 



Tau Lepton Production in ep Collisions at HERA 



HI Collaboration 



Abstract 



> 

(N 
(N 

O 

o 
\o 
o 

Q I The production of tau leptons in ep collisions is investigated using data recorded by the HI 

Q^' detector at HERA in the period 1994-2000. Tau leptons are identified by detecting their 

O . decay products, using leptonic and hadronic decay modes. The cross section for the pro- 

duction of tau lepton pairs is measured for the first time at HERA. Furthermore, a search for 
events with an energetic isolated tau lepton and with large missing transverse momentum is 
^ . performed. The results are found to be in agreement with the Standard Model predictions. 
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1 Introduction 

In the Standard Model (SM), tau leptons are produced either in pairs or in association with a tau 
anti-neutrino, as expected from lepton flavour conservation. In electron-proton collisions, pairs 
of tau leptons are produced via photon-photon interaction 77 -^ t^t^ (figure la), in which a 
photon from the electron interacts with a photon emitted by the proton [1]. Tau leptons and tau 
anti-neutrinos are produced in W boson decays, as illustrated in figure lb [2]. The signature of 
these events is a high transverse momentum (P^) isolated tau lepton, large missing transverse 
momentum P™'^^ due to the undetected neutrinos, and a hadronic system, typically of low P^. 

This paper presents a measurement of the production of tau lepton pairs {t^t^) and a search 
for events with an isolated tau lepton accompanied by large missing transverse momentum 
(r + P™'^^). The measurement of t^t~ production is performed at low transverse momentum 
considering both leptonic and hadronic tau decays. This measurement complements that of 
electron and muon pair production previously performed by the HI collaboration [3,4]. In 
the search for r + P™'^^ events tau decays are identified in the hadronic channel only. This 
search complements the previous HI measurements of the production of events with an isolated 
electron or muon and large missing transverse momentum, which have revealed an excess over 
the SM expectation of events containing in addition a high Pt hadronic system [5-7]. The 
ZEUS collaboration has also studied the production of events with an isolated lepton and large 
missing transverse momentum [8-10]. 

The analysis is based on data from electron ^-proton collisions at a centre-of-mass energy 
of 301 or 319 GeV, recorded by the HI experiment at HERA in the period 1994-2000. The 
total integrated luminosity amounts to 106 pb^^ for the measurement of t^t^ production and 
115 pb^^ for the search for r + P^^^^ events. 

This paper is organised as follows. The physics processes leading to tau lepton production at 
HERA are described in section 2 together with the relevant background processes. In section 3 
the HI detector and the experimental conditions are briefly described. Particle identification is 
presented in section 4. The selection of events with tau pairs and the resulting measurements 
are described in section 5. The search for events with an isolated high Pt tau lepton and large 
missing transverse momentum is presented in section 6. A summary is given in section 7. 



2 Signal and Background Processes 

The Standard Model processes leading to tau lepton production in electron-proton collisions 
are briefly outlined in this section, together with the dominant background processes. For each 
process, X is used to label the hadronic final state, excluding the tau decay products. 

The following processes, denoted henceforth as signal, lead to events containing genuine 
tau leptons in the final state: 



'in this paper, the name of the particle is used for both particles and anti-particles, unless otherwise stated (e.g. 
"electron" is used generically to refer to both electrons and positrons). 
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Figure 1 : Diagrams of the main production mechanisms of tau leptons in electron-proton col- 
lisions: a) tau pair production via photon-photon collisions and b) single W boson production 
followed by the subsequent decay of the W into a tau lepton and a tau anti-neutrino. 



• Tau pair production: ep -^ er^r^X 

The production of tau pairs proceeds mainly via photon-photon collisions [1], as shown 
in figure la. The proton can remain intact (elastic production, which dominates) or be 
dissociated in the interaction (inelastic production). The incident electron is usually scat- 
tered at small angles and is often not observed in the main detector. Only a small fraction 
of the total cross section is visible in the detector, as the cross section steeply falls with 
the transverse momentum P^ of the tau leptons. The cross section is about 20 pb for 
P^ > 2 GeV. This process is modelled using the GRAPE [11] Monte Carlo (MC) gener- 
ator. 

• Production of W bosons: ep -^ eWX -^ erz/^X 

The dominant production mechanism for W bosons in ep collisions [2] is depicted in 
figure lb. The cross section is largest in the photoproduction regime (photon-proton 
collisions), in which the W boson usually has small transverse momentum and the scat- 
tered electron is not observed in the main detector. The W production cross section 
at HERA is about 1 pb, resulting in a cross section of about 0.1 pb for the process 
ep — > eWX -^ erUrX. This process is modelled using the EPVEC generator [2]. The 
next-to-leading order QCD corrections to W production [12] are taken into account by 
weighting the events as a function of the rapidity and transverse momentum of the W 
boson [13]. 

The following processes, denoted henceforth as background, do not contain genuine tau 
leptons but contribute to the selected samples through misidentification or mismeasurement: 



Neutral Current deep inelastic scattering (NC DIS): ep -^ eX 

The scattered electron, or a quark or gluon that hadronises into a coUimated jet of low 
particle multiplicity, may fake the hadronic tau decay signature. Missing transverse mo- 
mentum may occur in NC DIS events because of fluctuations in the detector response and 



limited geometrical acceptance. The RAPGAP [14] generator is used to calculate this 
contribution to the background, including diffractive processes. 

• Photoproduction of jets: 7p ^ X 

Photoproduction processes may contribute to the background if a jet in the hadronic final 
state is misidentified as a hadronic tau decay. As in NC DIS events, missing transverse 
momentum may be measured due to fluctuations in the detector response and limited 
geometrical acceptance. The PYTHIA generator [15] is used to calculate the contribution 
from non-diffractive hard scattering photoproduction processes. The contribution from 
diffractive photoproduction processes is simulated using the RAPGAP generator. 

• Electron or muon pair production: ep —^ ee+e X or e/i+/i"X 

Events with two leptons (ee, fifi or e/i) may mimic events containing t^t^ pairs in which 
both tau leptons decay leptonically. The production of electron or muon pairs may also 
contribute as background if one electron or muon is not identified and fakes the signature 
of a hadronic tau decay. This contribution is calculated using the GRAPE generator. 

• Charged current deep inelastic scattering (CC DIS): ep —^ uX 

CC DIS events contain genuine large missing transverse momentum due to the presence 
of a neutrino in the final state. A jet originating from the struck quark may fake the 
hadronic tau decay signature. This process constitutes the main background contribution 
to the r+P™'^^ analysis. The DJANGO [16] generator is used to calculate the background 
contribution from CC DIS processes. 

All generated events are passed through the full GEANT [17] based simulation of the HI 
apparatus and are reconstructed using the same procedure that is applied to the data. 



3 Experimental Conditions 

At HERA, electrons of energy E^, = 27.6 GeV collide with protons of energy Ep = 820 or 
920 GeV, corresponding to a centre-of-mass energy ^/s of 301 or 319 GeV, respectively. The 
analysed datasets consist of 36.5 pb~^ of e~^p collisions at ^/s = 301 GeV (taken in the period 
1994-1997), 65.1 pb"^ of e+p collisions at v^ = 319 GeV (1999-2000) and 13.6 pb^^ of e'p 
collisions at ^/s = 319 GeV (1998-1999), corresponding to a total integrated luminosity of 
115.2 pb^^. The data recorded in the period 1994-1995 are not included for the measurement 
of tau pair production, reducing the integrated luminosity to 106.0 pb~^ for this analysis. 

A detailed description of the HI detector can be found in [18]. Only the components es- 
sential to the present analysis are described here. A right handed cartesian coordinate system 
is used with the origin at the nominal primary ep interaction vertex. The proton beam direction 
defines the z axis. The polar angles 9 and transverse momenta Pt of all particles are defined 
with respect to this axis. The azimuthal angle defines the particle direction in the transverse 
plane. The pseudorapidity is defined as ?] = — In tan |. 

The inner tracking system contains the central (20° < 9 < 160°) and forward (7° < ^ < 25°) 
drift chambers, and the proportional chambers which are employed for triggering purposes. It 



is used to determine the position of the interaction vertex and to measure the trajectories of 
charged particles. Particle transverse momenta and charges are determined from the curvature 
of the trajectories in a solenoidal magnetic field of 1.15 T. 

Hadronic final state particles as well as electrons and photons are absorbed in the highly 
segmented liquid argon (LAr) calorimeter [19] (4° < 9 < 154°) which is 20 to 30 radia- 
tion lengths deep depending on the polar angle. The hadronic section of the LAr calorimeter 
is 5 to 8 hadronic interaction lengths deep. Electromagnetic shower energies are measured 
with a precision of a{E)/E = 12%/ ^jE/GeV © 1% and hadronic shower energies with 
a{E)/E = 50%/^/E/GeV © 2%, as determined in test beam measurements [20]. In the back- 
ward region (153° < 9 < 178°), the LAr calorimeter is complemented^ by a lead-scintillating 
fibre spaghetti calorimeter (SpaCal) [21]. The LAr and SpaCal calorimeters are enclosed within 
the superconducting coil and are surrounded by an iron return yoke which is instrumented with 
streamer tubes to allow for the detection of muons in the range 4° < 9 < 171°. 

Dissociated proton states may be measured at small polar angles by a set of detectors in 
the forward direction: the PLUG, a sandwich calorimeter constructed from copper plates and 
silicon counters (0.6° < 9 < 3.5°), the proton remnant tagger (PRT), an array of scintilla- 
tors (0.06° < 9 < 0.17°), and the drift chambers of the forward muon detector (FMD) [23] 
(3° < ^ < 17°). 

The luminosity is determined from the rate of Bethe-Heitler events ep — > e'jp, where the 
photon is detected in a calorimeter located downstream of the interaction point. 

The events selected in this analysis are triggered by detecting electromagnetic clusters in the 
LAr or SpaCal calorimeter (electron trigger), by measuring a large missing transverse momen- 
tum in the LAr calorimeter (P^^^^^ trigger), or by using hits in the muon detectors combined with 
central tracker signals (muon trigger). In the kinematic range of this analysis, events contain- 
ing a pair of tau leptons are triggered with an efficiency of about 55% (35%) if one tau lepton 
decays leptonically (both tau leptons decay hadronically). Events containing single tau leptons 
and missing transverse momentum above 12 GeV (25 GeV) are triggered with an efficiency of 
about 50% (85%). 



4 Particle Identification 



4.1 Identification of electrons 

An electron candidate is defined by the presence of a compact and isolated electromagnetic 
cluster of energy deposits in the LAr or SpaCal calorimeter. The kinematics of the electron 
candidate are measured from the calorimeter cluster. Among the charged tracks reconstructed 
in the event, the track with the lowest extrapolated distance of closest approach to the cluster (6) 
is associated to the electron if it satisfies the condition 6 < 12 cm. In this case, the azimuthal 
angle and the charge of the electron are measured from the associated track. Additional energy 
within a cone of radius 0.5 in the pseudorapidity-azimuth {r]-(f)) plane around the electron 



In 1994, a lead-scintillator calorimeter [22] was used instead. 



candidate is required to be less than 3% of the energy attributed to the electron candidate. The 
electron identification efficiency is estimated using NC DIS events and is greater than 95% for 
an electron energy above 5 GeV. 



4.2 Identification of muons 

A muon candidate is identified by a track in the inner tracking system associated with a track 
segment or an energy deposit in the instrumented iron. The muon momentum is measured from 
the track curvature in the solenoidal magnetic field. The rejection of hadrons traversing the 
calorimeter and reaching the muon detectors is improved by requiring that the muon candidate 
deposits less than 5 GeV around its extrapolated track in the LAr calorimeter within a cylinder 
of radius 35 cm in the electromagnetic and 75 cm in the hadronic section. The efficiency to 
identify muons is estimated using elastic 77 -^ ji^jji^ events and is greater than 85% in the 
energy range considered in this analysis. 



4.3 Reconstruction of the hadronic final state 

The hadronic final state (HFS) is measured [24] by combining energy measurements from the 
calorimeter with charged particle momenta measured by the inner tracker. Identified isolated 
electrons or muons are excluded from the HFS. The hadronic energy scale is calibrated by 
comparing the transverse momentum of well measured electrons to that of the HFS in NC DIS 
events [25]. Jets with P-r > 2 GeV are reconstructed from HFS particles using an inclusive kx 
algorithm [26] in the Pt recombination scheme with a separation parameter set to one. 



4.4 Identification of tau decays 

Leptonic tau decays (branching ratio 35% [27]) are identified by detecting an electron or a 
muon as described above. Hadronic tau decays typically produce low multiplicity, coUimated, 
hadronic jets, henceforth denoted as r-jets. Depending on the number of charged hadrons 
produced, the hadronic decay modes are summed up in two categories, referred to as "1-prong" 
(one charged hadron, branching ratio 49%) and "3-prong" (three charged hadrons, branching 
ratio 14%). The branching ratio for decays into more than three charged hadrons is small (about 
2%) and such decays are not considered in the present analysis. 

Two different algorithms to identify hadronic tau decays, applied to jets reconstructed in the 
angular range 20° < 9 < 120°, have been developed and are described below. The measurement 
of tau lepton pair production, in which the tau leptons generally have low momentum, requires 
an optimal background rejection. In contrast, the search for tau leptons produced in W decays 
at high Pt uses an algorithm that maximises the identification efficiency, since the signal cross 
section is low and the background is less severe. 



A: Neural Network based identification algorithm optimised for low Pt taus 

For the measurement of the tau lepton pair production, an algorithm has been developed that 
uses multiple neural networks to discriminate between r-jets and the background from elec- 
trons, muons or hadronic jets. A detailed description of the algorithm can be found in [28]. The 
algorithm is implemented in two steps. 

In the first step, 1-prong (3-prong) candidates for hadronic tau decays are preselected by 
requiring exactly one (three) well reconstructed track(s) in the drift chambers within a cone 
("r-cone") around the jet axis. The opening angle of the r-cone varies between 5° and 30° 
depending on the jet momentum, with smaller angles at higher momentum due to the larger 
Lorentz boost in the direction of the tau candidate. The tracks are required to be not associated 
to identified electrons and muons and the scalar sum of their transverse momenta is required to 
be larger than 2 GeV. The fine granularity of the LAr calorimeter is used to match extrapolated 
tracks with energy deposits in the calorimeter and to reconstruct additional neutral particles 
associated to the tau candidate from unmatched energy deposits in the r-cone. The sum of 
the four-vectors of the tracks and of the neutral particles defines the r-jet four-vector. If all 
associated tracks have a well measured charge, the charge of the tau candidate is reconstructed 
as the sum of the charges of the associated tracks. 

In the second step of the algorithm, various variables related to the particle multiplicity and 
coUimation of the r-jet candidate are used. The set of variables includes: the multiplicity of the 
neutral clusters within the r-cone; the invariant masses calculated from clusters, from tracks 
and from charged and neutral particles in the r-cone; the difference in energy measured from 
tracks and from clusters; the distance ini] — (p between tracks and clusters; the first and second 
transverse moments of the distribution of energy deposits in the calorimeter with respect to the 
jet axis; the sum of energy deposits detected in an extended cone of radius 1.0 in ?7-0 around 
the r-cone. Two neural networks (NN) using these variables are separately trained using MC 
simulations to identify 1-prong and 3-prong tau decay modes, the output of which is denoted 
by /,1-prong and /,3-prong^ respectively. Their numerical value varies between zero and one and 
is used to discriminate between tau candidates (close to one) and hadronic jets (close to zero). 
Depending on the track multiplicity of the jet, the output of one of these NNs is used to select 
tau candidates. By requiring L^-p^°^s > o.75 (for 1-prong candidates) or /,3-prong ^ 0^5 (foj- 
3-prong candidates), the efficiency of this algorithm to identify hadronic tau lepton decays in 
T^T^ events is about 50%, as calculated from MC simulations of tau decays with visible energy 
in the range considered in this analysis. The probability for hadronic jets to be misidentified as 
hadronic 1-prong (3-prong) tau decays is 0.5% (4%). 

The signature of hadronic 1-prong tau decays may also be faked by unidentified electrons 
and muons. Two additional neural networks are trained to veto such cases. The output of 
these NNs, denoted by -^^cto ^^^ -^vcto respectively, is expected to be close to one for tau de- 
cays and close to zero for electrons or muons that fail the identification criteria described in 
sections 4.1 and 4.2. 

B: Identification algorithm optimised for high Pt taus 

A different approach to tau identification is used to search for events with a high momentum 
tau lepton and large P™'^^. A high identification efficiency and a sufficient level of background 



rejection are achieved by requiring a collimated jet, containing only one charged particle and 
isolated with respect to other tracks and jets within a cone of radius 1.0 in rj-cj). The identifica- 
tion of hadronic r-decays is therefore restricted to hadronic 1-prong decay modes. 

The coUimation of a jet is measured by the jet radius Rjet, defined as the energy weighted 
average distance in 7^-0 between the jet axis and all HFS particles h (neutral and charged) in 
the jet: 

^^•-* = ^y2 EhV^vijet, hf + A0(jet, hf . (1) 

Ejet ^ 

Jets with one track and Rjet < 0.12 are selected as tau candidates. The four-vector of the tau 
candidate is taken to be that of the jet. For jets with Pt > 7 GeV, this identification procedure 
results in an efficiency of about 80% to identify 1-prong hadronic decays of tau leptons resulting 
from decays of W bosons. The misidentification probability for hadronic jets is less than 1%. 



5 Production of r+r Pairs 



5.1 Event selection 

The decay modes investigated in the present analysis are classified as leptonic when both taus 
decay leptonically (branching ratio 6.2%, excluding decays to same flavour leptons), 
semi-leptonic when one tau decays hadronically and the second leptonically (45%) and hadronic 
when both taus decay hadronically (42%). The case where both tau leptons decay to charged 
leptons of the same flavour (ee or jiji final state) is not included in the present analysis, as the 
separation of the t^t^ signal from electron or muon pair production is difficult. 

Leptonic and hadronic tau decays are identified as described in the previous section and by 
applying in addition the following selection criteria. The isolation of electrons and muons is 
measured by the distance in the r^ — plane to the closest hadronic jet (-Djet) and to the closest 

track (Dtrack)- 

• e candidates are reconstructed in the polar angular region 20° < O^, < 140° and are 
required to have an energy above 5 GeV and a transverse momentum above 3 GeV. They 
must be isolated from jets with D^^t > 10. 

• /i candidates are reconstructed in the polar angular region 20° < 9^ < 140° and are 
required to have a transverse momentum above 2 GeV. They must be isolated from other 
tracks and jets with -Dtrack > 0.5 and Djct > 1-0, respectively. 

• r-jet candidates of transverse momenta above 2 GeV are reconstructed in the polar an- 
gular region 20° < 9r < 120°, using algorithm A as described in section 4.4. The output 
of the neural network is required to satisfy L^-p^°''^ > 0.75 or L^"P''°°s > 0.75. 
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Events with two tau candidates are selected. If the charge of both tau candidates is measured, 
events with two candidates of the same charge are rejected. 

In order to avoid significant background contributions from NC DIS and 7p processes, the 
analysis is restricted to elastic ep -^ epr^r^ production. Inelastic events are vetoed by re- 
quiring no extra track or energy deposit above the noise level in the main detector in addition 
to those associated to the decay products of both tau leptons and a possible additional elec- 
tron. Furthermore, no significant activity should be observed in the forward detectors (PLUG, 
FMD and PRT). The requirements applied [29] ensure that the proton remains intact or disso- 
ciates into a low mass state. Remaining background originates mainly from elastic lepton pair 
production and from diffractive NC DIS or -jp processes. 

Electron and muon pair production processes (77 — * e~^e~ and 77 -^ /U+yU") constitute a 
background to tau pair production in semi-leptonic decay modes when one of the leptons is 
correctly identified as e or yU while the second lepton fails to be identified by the algorithms 
described in sections 4.1 and 4.2 and fakes the signature of a 1-prong hadronic tau decay. In 
order to reject this background, L^^^^ > 0.75 (L^cto > 0.75) is required for 1-prong tau jet 
candidates selected in the e + r-jet (/i + r-jet) channel. 

To further reduce the background from NC DIS processes, where the scattered electron is 
selected as a candidate for a r ^ e decay and the struck quark fragments into a coUimated jet 
of low particle multiplicity that fakes the r-jet signature, the longitudinal momentum balance 
calculated as 

E-P, = ^E,(l-cos^,) (2) 

i 

is employed, where Ei and 9i denote the energy and polar angle of each particle detected in the 
event. For events in which only momentum in the proton direction is undetected, E — P^is equal 
to twice the energy E^. of the electron beam, i.e. 55.2 GeV. For events containing tau leptons, 
E — Pz values well below 55.2 GeV are expected since the neutrinos produced in the tau decays 
are not detected and, most of the time, the scattered beam electron escapes down the beam pipe. 
By requiring _E — P^ < 50 GeV when the detected electron has the same charge as the beam 
lepton, the NC DIS background is rejected to a large extent. The E — P^ cut is also applied to 
the leptonic channel and rejects muon pair-production events, for which the scattered electron 
is detected together with one produced muon while the second muon escapes in the forward 
direction. If a second electron is detected in the event, this condition is not applied. 

The selection criteria for t^t^ events are summarised in table 1. With these selection 
criteria, 1.2% of the elastic 77 -^ t^t^ events in which both tau leptons satisfy Pf > 2 GeV 
and 20° < 9^ < 140° are selected. The efficiency is limited by the fact that the energy of the 
detected tau decay products is significantly lower than the tau lepton energy, since the neutrinos 
escape detection. 



5.2 Background studies 

Due to the aforementioned elasticity requirements, the remaining background from NC DIS 
and 7p processes consists mainly of exclusive diffractive events, for which the validity of the 
resolved pomeron model [30] as implemented in the RAPGAP program is questionable. Hence 
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Selection of r+r Events 


Decay Channel 


Leptonic 


Semi-leptonic 


Hadronic 


r Signatures 


e 


e 
r-jet 


r-jet 


r-jet 
r-jet 


Events with two same charge tau candidates rejected 


Elastic 
Production 


No additional tracks, no additional clusters in the LAr or SpaCal calorimeter, 
no activity in forward detectors above noise level 


Background 
Reduction 




L^eto > 0.75(1) 


^v^cto > 0.75(1) 




E-P,< 50 GeV(2) 






(^) applied only to 1-prong r-jet candidates 

(^) applied only if the electron associated with the tau decay has the same charge as the beam lepton 
and no second electron is detected 



Table 1 : Selection criteria for elastic t^t events in the leptonic (e + /i), semi-leptonic (e + r- 
jet, /i + r-jet) and hadronic (r-jet + r-jet) decay modes of the r lepton pair. 



the RAPGAP prediction for NC DIS (7p) is normalised to the number of events observed in a 
control sample in which an electron and a jet (two jets) are selected in the Pt and 9 ranges of 
the analysis, and where the elastic requirements are applied. It has been verified that the shapes 
of the observed kinematic distributions are reasonably well described by RAPGAP. 

Furthermore, to check that the probability for an electron, a muon or a hadronic jet to be 
misidentified as a r-jet candidate is well described by the MC simulation, several event samples 
are studied in which the contribution of individual background processes is enhanced. These 
event samples are selected in a phase space similar to that of the t^t^ event sample, requiring 
the conditions for elastic production to be fuUfilled. In order to enhance the background com- 
ponent in the control samples, the condition L^-p^°''^ > 0.75 or L^"P''°°s > 0.75 on the NN 
outputs for r-jet candidates is not applied. 

The following samples are employed: 



77 -^ e+e~ control sample 

A 77 -^ e+e^ dominated event sample is defined by selecting events with one electron 
and one 1-prong r-jet candidate. No veto condition on L^cto i^ applied. Events with two 
same charge tau candidates are rejected. 

77 -^ At^A*" control sample 

A 77 -^ ^+^- enriched event sample is defined by selecting events with one muon and 
one 1-prong r-jet candidate. No veto condition on Lt^cto i^ applied. Events with two 
same charge tau candidates are rejected. 
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• 7P control sample 

A 7p enriched event sample is defined by selecting events containing two jets with one, 
two or three tracks. If both jets are 1-prong tau candidates, the requirements L^^^^ > 0.75 
and L^eto > 0-75 are applied to suppress the contributions from 77 -^ e+e~ and 77 -^ 
fi~^H~ processes. In order to reduce the contribution from NC DIS process, only events 
with E ~ Pz < 4:5 GeV are accepted. 

• NC DIS control sample 

A NC DIS enriched event sample is defined by selecting events with one electron and one 
r-jet candidate. If the r-jet is a 1-prong candidate, L^cto > O'^'S is required in order to 
suppress the contribution from 77 -^ e+e^ processes. 

The selection criteria of the background control samples are summarised in table 2, in which 
the numbers of events obtained from the data and the MC simulation are also given. 

The distributions of all quantities used in the selection of the t^t^ event sample are well 
described by the MC simulation both in shape and normalisation in the control samples. Ex- 
amples of these distributions are shown in figure 2. The distributions of the electron and muon 
rejection discriminators L^^to ^^^ -^vcto ^^ shown in figures 2a and 2b for the control samples 
where most of the tau candidates are unidentified electrons and muons, respectively. The distri- 
bution of the NN output L^^p^^s (^s-prong^ |g gj^Q^j^ jj^ fjg^j.g 2c (2d) for 1-prong (3-prong) 
jets in the ^p control sample which is enriched in hadronic jets. In figures 2a-d the contribution 
from tau pair production populates the region close to one, while the background accumulates 
at values close to zero, as expected. Finally, the distribution of the E — P^ variable is shown 
in figure 2e for the NC DIS control sample. The agreement between data and simulation in the 
control samples shows that the background contribution as well as the experimental efficiencies 
are modelled by the MC simulation within the attributed systematic uncertainties described in 
the next section. 



5.3 Systematic uncertainties 

In this section, the systematic uncertainties associated with the measurement of elastic t^t^ 
production are discussed. The uncertainties on the signal expectation are determined by varying 
the following experimental quantities by ±1 standard deviation. 

• Electron identification and reconstruction 

In the kinematic region considered in the t^t^ analysis, the uncertainty on the electron 
identification efficiency is 5% (2% when the electron energy is above 10 GeV). The elec- 
tron energy scale uncertainty is estimated to be 3%. The uncertainties on the measurement 
of the electron angles 9 and (p are 3 mrad and 1 mrad, respectively. 

• Muon identification and reconstruction 

The uncertainty on the muon identification efficiency is 5%. The relative scale uncer- 
tainty on the muon momentum is conservatively taken to be 5%. The uncertainties on the 
measurements of the muon angles 9 and are 3 mrad and 1 mrad, respectively. 
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T~^T Background Control Samples 


Control 
samples 


e+e 


fi+fi 


IP 


NCDIS 


Signatures 


e 
1-prong r-jet 


1-prong r-jet 


two jets 
with < 3 tracks 


e 
r-jet 


Events with two same charge 
tau candidates rejected 






Elastic 
Production 


No additional tracks, no additional clusters in the LAr or SpaCal calorimeter, 
no activity in forward detectors above noise level 


Background 
Reduction 






L^eto > 0.75(1) 

^Ceto > 0.75(1) 

E-P,<A5 GeV 


L^^to > 0.75(2) 


HI Data 
Total SM 


115 
133.1 ±19.5 
(95% e+e") 


20 

14.1 ±1.4 

(50% n+fi-) 


29 

24.9 ± 10.0 

(79% 7p) 


29 

32.4 ±6.3 
(62% NC DIS) 


(^) applied only if both jets are 1-prong tau candidates 
(^) applied only to 1-prong r-jet candidates 



Table 2: Selection criteria for the background control samples in which each background contri- 
bution is individually enhanced. The numbers of observed and expected events are also shown. 
The dominant contribution to the total SM expectation is indicated as a percentage in the last 
row. Here, the tau candidates are not required to satisfy the condition L^~p^°^s > 0.75 or 

2^3-prong > 0.75. 



Identification and reconstruction of hadronic tau decays 

For each of the charged tracks associated to the r-jet, a reconstruction efficiency un- 
certainty of 3% is assigned. The efficiency to identify hadronic r decays with the NN 
algorithm A has an additional uncertainty of 10% [28], estimated by comparing different 
simulations of shower development in the LAr calorimeter. It has been verified, using 
a large statistics sample of hadronic jets from inclusive CC DIS and NC DIS samples, 
that the output of both NNs is well described by the MC simulation. The uncertainties 
arising from the modelling of tau decays are negligible, as estimated by comparing the 
results obtained using either the PYTHIA [3 1] or TAUOLA [32] programs to simulate the 
decays of the tau leptons generated with GRAPE. The energy of the neutral clusters of 
r-jet candidates has a relative uncertainty of 4%, corresponding to the hadronic energy 
scale uncertainty. The uncertainties on the measurements of the r-jet angles 9 and are 
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estimated to be 10 mrad and 5 mrad, respectively. 

• Elastic event selection 

The efficiency with which elastic events are selected and proton dissociative events are 
rejected depends on the noise level in the LAr and on the performance of the forward 
detectors [29]. Its uncertainty does not exceed 3%. 

• Triggering 

The trigger efficiency is studied using elastic e~^e~ and fi'^fi~ events and diffractive -jp 
di-jet events in a phase space similar to that of the t^t^ analysis. The uncertainty on 
the trigger efficiency depends on the region in 9 in which the tau candidates are detected: 
central (9 > 30°) or forward (9 < 30°). The uncertainty is 10% if both tau candidates are 
in the central region, 20% if one tau candidate is detected in the forward and the other in 
the central region and 30% if both tau candidates are detected in the forward region. 

• Luminosity 

The luminosity of the analysed datasets is measured with an uncertainty of 1.5%. 

The individual effects of the above experimental uncertainties are combined in quadrature, 
yielding a total uncertainty of 21% on the signal expectation. The largest contributions to this 
uncertainty arise from systematics attributed to the tau identification procedure (15%) and to 
the trigger efficiency (12%). 

Contributions from background processes, modelled using the generators described in sec- 
tion 2, are attributed relative systematic uncertainties of 50% (7p), 30% (NC DIS), 
15% (77 -^ e+e~, 77 -^ fi^n^), estimated from the level of agreement observed between 
the MC simulation and the data in the background enhanced control samples described in sec- 
tion 5.2. 



5.4 Results 

In total, 30 t^t^ candidate events are selected, in agreement with the SM expectation of 
27.1 ± 4.1 events, of which 16.0 ± 3.4 are expected from 77 -^ t^t~ signal processes. The 
signal expectation is dominated (85%) by the elastic production component. The lepton pair 
production, NC DIS and -jp processes contribute with similar rates to the background expecta- 
tion. The numbers of observed and expected events in the four analysed channels are shown in 
table 3. An event selected in the semi-leptonic channel is displayed in figure 3. 

The distributions of the polar angle and of the transverse momentum of the tau candidates, 
together with the visible invariant mass, are shown in figure 4. The measured distributions are 
compatible with the SM expectations. As expected for 77 -^ t^t^ processes, most tau decay 
products are detected with relatively small transverse momenta. 

Using the selected sample, a measurement of the cross section for the elastic production 
of T~^T^ pairs is performed in the kinematic region defined by 20° < 9^- < 140° and P^ > 
2 GeV. For this measurement, the data samples collected at ^/s = 301 GeV and 319 GeV are 
combined, taking into account their respective luminosities. Assuming a linear dependence of 
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T~^T Results 


Decay Channel 


Leptonic 


Semi-leptonic 


Hadronic 

r-jet r-jet 


Total 


e r-jet 


fj, r-jet 


HI Data 
SM 

r+r- 


^7 

1 

2.9 ±0.4 

56% 


2 
6.3 ±0.9 

47% 


10 
7.0 ±1.3 

85% 


11 

11.0 ±2.0 

50% 


30 

27.1 ±4.1 
59% 



Table 3: Number of selected events and SM prediction for the t^t analysis. The expected 
relative contribution of the r^r^ process to the SM prediction is also shown. 



the cross section on the proton beam energy, as predicted by the SM, the measured cross section 
corresponds to an effective centre-of-mass energy of ^/s = 314 GeV. The cross section is 
calculated using the formula: 

-''data -''bgr ^on 

" = C-A ' ^'^ 

where A^data is the number of observed events, A^bgr the expected contribution from background 
processes, C the total integrated luminosity and A the signal acceptance. The contribution 
from inelastic 77 — > t'^t^ processes is included in the background expectation. The signal 
acceptance A is calculated using the GRAPE generator, as the ratio of the number of events 
accepted at reconstructed level to the number of events generated in the defined phase space. It 
accounts for the selection and trigger efficiencies and for the differences in momentum between 
the original r leptons and the detected decay products. 

The measured cross section for elastic tau pair production ep — * epr^r^ integrated over 
the phase space defined above is 13.6 ± 4.4 ± 3.7 pb where the first error is statistical and 
the second systematic. The result is in agreement with the SM expectation of 11.2 ± 0.3 pb, 
calculated using the GRAPE generator. 



6 Production of High Pt Tau Leptons in Events with Large 
Missing Transverse Momentum 

6.1 Event selection 



Events containing an isolated tau lepton and large missing transverse momentum are selected 
with a procedure similar to that used in the analysis of events with an isolated electron or muon 
and large missing transverse momentum [7]. The tau leptons are identified using hadronic 
decays only, as the leptonic tau decays lead to final states which cannot be distinguished from 
those studied in [7]. 

The event selection is performed in two steps. In the first step, the net transverse momentum 
reconstructed from all particles (electrons, muons and hadrons) P™'^^, is required to be above 
12 GeV. In order to ensure uniform trigger conditions, the net transverse momentum measured 
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from all energy deposits detected in the calorimeter, P^^^°, is required to be above 12 GeV. 
The reconstructed P™'^^ is approximately equal to P|.^'° except for events containing muons 
in the final state. In order to exploit further the event topology in the transverse plane, the 
variable Vap/Vp is employed, defined as the ratio of the anti-parallel to parallel projections of 
all energy deposits in the calorimeter with respect to the direction of P^^^° [33]. Events with 
genuine missing transverse momentum are in general reconstructed with Vap/Vp values close 
to zero and large values of P|.^'°, whereas background events from NC DIS and -jp processes 
are intrinsically balanced, with larger values of Vap/Vp and low values of P^^^°. Only events 
with Vap/Vp below 0.5 are accepted. For events with P|,^'° below 25 GeV, a stricter condition 
of Vap/Vp < 0.15 is applied. The background from NC DIS is further reduced by requiring 
_E — P^ < 50 GeV. Additionally, the events are required to contain at least one isolated hadronic 
jet with transverse momentum above 7 GeV in the central region of the detector 
20° < 9 jet < 120°. The isolation is characterised by the distance of the jet in r^ — to the 
nearest hadronic jet (-Djct > 10) and the nearest track not belonging to the jet (-Dtrack > 10). 

In the second selection step, the isolated jets are required to pass the tau identification criteria 
of algorithm B described in section 4.4. A jet is accepted if it is narrow (Rjet < 0.12) and 
contains exactly one charged track (A^/^acfcs — !)• ^^ order to remove background from hadronic 
jets containing mostly neutral particles, this track is required to have a transverse momentum 
ptrack gj-eater than 5 GeV. If more than one isolated jet satisfies these requirements, the one with 
the highest Ft is considered as the tau candidate. In order to further reduce the background from 
intrinsically balanced events, in which the P™'*^*^ is due to a mismeasurement, the acoplanarity 
A0, defined as the angle in the transverse plane between the r-jet candidate and the hadronic 
system excluding the tau candidate (X), is required to be below 170°. This criterium removes 
events in which the r-jet candidate and the rest of the hadronic system are back-to-back, as is 
typical for NC DIS and 'jp events. 

A summary of all selection criteria is presented in table 4. Using these selection criteria, 
SM W ^ Tu events are selected with an overall efficiency of 8%. In comparison to W decays 
into electrons or muons [7], 14^ ^ ru decays are selected with a significantly lower efficiency, 
mainly due to the branching ratio (49%) for hadronic 1-prong tau decays and to the more 
restricted polar angular range of this analysis. 



6.2 Background studies 

After applying the selection criteria, the main background is expected to occur from events with 
genuine missing transverse momentum produced by CC DIS processes, in which a narrow jet 
with low track multiplicity fakes the tau signature. Additional background arises from NC DIS 
and 7p processes, which have a much larger cross section than the CC DIS process and lead to 
events that contain hadronic jets but no genuine P™''^^ The modelling of the CC DIS, NC DIS 
and 7p backgrounds is verified using samples in which the contribution of each background 
process is enhanced. The selection criteria defining the background control samples are listed 
in table 5, together with the observed number of events and the SM expectation. 

The CC DIS background control sample is selected using the "P^'^s _,_ isolated jet" prese- 
lection described in section 6.1. The NC DIS and 7p samples are selected in a complementary 
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Selection of r + P™'^^ Events 




ocalo 


> 


12GeV 




pmiss + 

Isolated Jet 


pmiss 


> 


12GeV 




Preselection 


E-P, 


< 


50GeV 






VaJVp 


< 


0.5 (< 0.15 for P|^'° 


< 25 GeV) 




Njets 


> 


1 






pjet 


> 


7GeV 






20° 


< 


Ojet < 120° 






^track 


> 


1.0 






Djet 


> 


1.0 






tracks 


= 


1 




Final 


JDtrack 


> 


5GeV 




T + P^'^" 


Rjet 


< 


0.12 




Selection 


A0 


< 


170° 





Table 4: Selection criteria for events containing an isolated r lepton and large P™' 



phase space at large E — P^ and large Vap/Vp, respectively. In these background samples the 
tight selection criteria of tau candidates are not applied and only isolated jets, as defined in 
table 4, are considered. 

As shown in figure 5a-c, the distributions of Rjet, ^l^^^^^, and P^'^'^^ in the CC DIS control 
sample are well described by the MC simulation, both in shape and normalisation. The dis- 
tribution of the acoplanarity angle A0 in the control sample enriched in 7p events, shown in 
figure 5d, shows a clear peak towards 180°, corresponding to back-to-back events. In figure 5e 
the distribution of the hadronic transverse momentum P^ for the NC DIS control sample is 
shown. The good agreement between data and MC simulation observed in all control samples 
confirms the good understanding of background rates and of the properties of the jets (shape 
and multiplicity) used in the tau identification procedure. 



6.3 Systematic uncertainties 

In this section, the systematic uncertainties associated with the search for events containing an 
isolated tau lepton and large missing transverse momentum are discussed. The effect of those 
uncertainties on the expectations from SM W production and from background processes is 
determined by varying the experimental quantities by ±1 standard deviation. 
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T + P^^^^ Background Control Samples 


Control 
Sample 


CCDIS 


NCDIS 


IP 


Selection 


P^'^^ > 12 GeV, P|^'° > 12 GeV, at least one isolated jet 


E-P,< 50 GeV 

Vap/Vp < 0.5 

Vap/Vp < 0.15 for P|^i° < 25 GeV 


E - P^ > 35 GeV 
P| > 10 GeV 

5° < ^^ < 150° 


E-P,<50 GeV 
VjVp > 0.15 
no electrons 


Data 
Total SM 


1811 

1858.6 ±120.7 

(93% CC) 


108 

100.3 ±22.0 

(98% NC) 


1165 

1173.1 ±208.3 
(80% 7P) 



Table 5: Selection criteria for the r ± P^^^^ background control samples together with the 
number of selected events compared to the SM expectation. The dominant contribution to the 
SM prediction is indicated as a percentage in the last row of the table. 



Tau identification procedure 

The main experimental uncertainty on the signal and background expectations arises from 
the requirement of exactly one track within the jet and from the condition on the jet radius 
Rjet < 0.12. The energy of charged tracks is reconstructed with an uncertainty of 5%. 
The measurements of the polar and azimuthal angles are attributed a systematic error of 
3 mrad and 1 mrad respectively. The uncertainty on the track reconstruction efficiency 
is 3% and has a relatively large effect on the background expectation, due to migrations 
of hadronic jets with multiple charged particles into the single track category. The mod- 
elling of Rjet is studied with high statistical precision using jets in a dedicated inclusive 
NC DIS sample. The measured and simulated jet radius distributions are compared and 
an uncertainty depending on the polar angle of the r-jet is attributed to Rjet- The uncer- 
tainties associated with the tau identification are 15% for the expected W signal and 25% 
for the background predictions. 

Hadronic final state reconstruction 

The hadronic energy measurement has a relative uncertainty of 4%. The polar angle 
measurement has an uncertainty varying from 3 mrad for jets reconstructed in the central 
region to 10 mrad for forward jets. The topological variable Vap/Vp is attributed a relative 
uncertainty of 10%. 

Trigger 

The uncertainty on the P^^^^ trigger efficiency is 5% deduced from a large statistics NC 
DIS sample for which the trigger information is reconstructed offline, ignoring the signal 
from the electron [33]. 
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• Luminosity 

The luminosity of the analysed datasets is measured with an uncertainty of 1.5%. 

• Theoretical errors on signal and background contributions 

The uncertainty on the W production signal cross section is estimated to be 15% [12]. 
The expectations from NC DIS and ^p processes are each attributed an additional uncer- 
tainty of 20%, a value which has been estimated from the control samples described in 
section 6.2. This uncertainty covers the sensitivity of the jet radius and multiplicity (used 
in the tau identification algorithm) to the modelling of parton showers in NC DIS and 7p 
MC samples. 

The individual effects of the experimental uncertainties are combined in quadrature to give 
the total experimental systematic uncertainty. The total uncertainty on the SM W signal is 20% 
and that on the SM background is 34%. For both signal and background processes, the total 
uncertainty is dominated by the uncertainty arising from the tau identification procedure and 
the theoretical uncertainties. 



6.4 Results 

In the final event sample 6 events are observed in the data, compared to a total SM expectation 
of 9.9 tl'^ events, of which 0.89^q 26 are expected from SM W production. The P^ spectrum 
and other properties of the events in the final sample are shown in figure 6. Table 6 summarises 
the results. The events observed in the data are concentrated in the region of very low P^ , 
where the contribution from CC DIS background processes dominates the SM expectation. In 
the region P^ > 25 GeV, where an excess of events containing isolated electrons or muons is 
observed [7], no event is found for a SM prediction of 0.39 ± 0.10, of which 0.20 ± 0.04 are 
expected from SM W production. 

In the absence of a signal, a model independent upper limit on the cross section for the 
production of events containing an isolated tau leptons and large missing transverse momentum 
is set in the kinematic region: 5° < ^^ < 140°, Pf > 10 GeV and P™'"" > 12 GeV. The 
limit is calculated using a modified frequentist approach based on likelihood ratios [34] and 
taking into account the systematic uncertainties discussed in section 6.3. The acceptance for 
processes producing isolated tau leptons in events with large missing transverse momentum 
in the given kinematic region is estimated using the MC simulation for SM W production, 
implemented in the EPVEC generator. An additional model uncertainty of 10% is attributed to 
the acceptance. This uncertainty is estimated by comparing the acceptance predicted by EPVEC 
with that obtained using the generator ANOTOP, which simulates the anomalous production of 
single top quarks in ep collisions at HERA [35] and produces W bosons with different kinematic 
distributions^. 

An upper limit of o- < 0.85 pb at 95% confidence level is obtained for the production 
cross section of events containing an isolated tau lepton and large missing transverse mo- 
mentum in the phase space defined above. In the region P^ > 25 GeV, the upper limit is 

^The P^ distribution for top events is peaked at about 70 GeV, while events from SM W production are 
concentrated at low P^ values. 
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T + P^'"" Results 


HI Data 


SM Expectation 


SM Signal 


Other SM Processes 


Total 


6 


q q +2.5 


0.89 t'oil 


9.0 til 


P# > 25 GeV 





0.39 +°:?? 


0.20 in^ 


0.19 IH^ 


Pf > 40 GeV 





n 1 fi +o-o'^ 

U.iD _Q06 


0.08 in^ 


0.08 l°:°e^ 



Table 6: Number of events observed in the data compared to the SM expectation for signal (W 
decay into r) and background processes in the final r + P™'^^ sample. 

(t(P^ > 25 GeV) < 0.31 pb at 95% confidence level. These limits are higher than the cross 
sections measured in the electron and muon channels [7]. The present measurement is therefore 
compatible with the previous measurement of events with an electron or muon and P^^^^, as 
expected if lepton universality is assumed. 



7 Summary 

In this paper, the production of tau leptons in ep collisions at HERA is investigated in events 
containing a r^r^ pair and events containing an isolated tau lepton and large missing transverse 
momentum. 

The production cross section of t^t^ pairs is measured in elastic processes, in a combi- 
nation of leptonic, semi-leptonic and hadronic decay modes of the two tau leptons. In a data 
sample corresponding to an integrated luminosity of 106 pb^\ 30 events are observed, in agree- 
ment with a Standard Model expectation of 27.1 ± 4.1 events, of which 16.0 ± 3.4 are expected 
from 77 -^ t^t^ signal processes. This is the first observation of tau pair production in ep 
collisions. 

A search for the production of isolated tau leptons in events with large missing trans- 
verse momentum is performed in a data sample corresponding to an integrated luminosity of 
115 pb^^. The selection yields 6 candidate events, compatible with a Standard Model expecta- 
tion of 9.9 I3 g events. No event is observed in the region P^ > 25 GeV, where 0.39 ± 0.10 
events are expected, including 0.20 ± 0.05 events from W -^ tv^. decays. An upper limit on 
the production cross section of cr[P^ > 25 GeV) < 0.31 pb at 95% confidence level is set. 
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Figure 2: Distributions of a) the likelihood Lycto for ^-jet candidates in the 77 -^ e^e^ control 
sample, b) the likelihood L'^^^^ for r-jet candidates in the 77 -^ /U^/U^ control sample, c) 
and d) the NN outputs L^~v™^^ and L^-prong ^^^ i-prong and 3-prong r-jet candidates in the 
7P control sample, and e) the longitudinal momentum balance E — P^ in the NC DIS control 
sample. In each figure the open histogram shows the total SM expectation and the shaded band 
its uncertainty. The contribution from tau pair production is shown as the hatched histogram. 
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Figure 3: A tau pair candidate event with one tau lepton decaying leptonically to a muon and 
the other tau lepton decaying to three charged hadrons (3-prong topology). 
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Figure 4: Distributions of a) the polar angle and b) the transverse momentum of the recon- 
structed r lepton candidates in the t'^t^ event sample. The invariant mass of the tau candidate 
pair is shown in c). The selected events enter distributions a) and b) twice. In each figure 
the open histogram shows the total SM expectation and the shaded band its uncertainty. The 
contribution from tau pair production is shown as the hatched histogram. 
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Figure 5: Control distributions in the background enriched samples defined for the r + P™'^^ 
analysis: a) the radius of isolated jets, b) the number of charged particles in each isolated jet and 
c) the transverse momentum of the track associated to single track isolated jets in the CC DIS 
control sample; d) the acoplanarity angle between the isolated jet and the remaining hadronic 
system X in the 7p control sample; e) the transverse momentum of the hadronic system X 
excluding the isolated jet candidate in the NC DIS control sample. In all distributions each 
event may enter several times, depending on the number of isolated jets in the event. In each 
figure the data (points) are compared to the SM expectation (solid histogram) shown with its 
uncertainty (shaded band). The signal contribution dominated by the SM W production is also 
shown (dashed histogram). In a) and b) it is shown multiplied by a factor 100 and in d) by 200. 
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Figure 6: Distributions of a) the total missing transverse momentum P^^^^, b) the hadronic 
transverse momentum P^ , c) the polar angle 9^. and d) the transverse momentum Pf of the 
r-jet candidate in the final sample of events containing an isolated tau lepton and large missing 
transverse momentum. In each figure the data (points) are compared to the SM expectation 
(solid histogram) shown with its uncertainty (shaded band). The signal contribution dominated 
by the SM W production is also shown (dashed histogram). 
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